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Abstract: 

We construct the Baxter Q— operator and the representation of the Separated Variables (SoV) 
for the homogeneous open SX(2,R) spin chain. Applying the diagrammatical approach, we 
calculate Sklyanin's integration measure in the separated variables and obtain the solution to 
the spectral problem for the model in terms of the eigenvalues of the Q— operator. We show that 
the transition kernel to the SoV representation is factorized into the product of certain operators 
each depending on a single separated variable. As a consequence, it has a universal pyramid-like 
form that has been already observed for various quantum integrable models such as periodic 
Toda chain, closed SL(2, R) and SL(2, C) spin chains. 
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1. Introduction 

Recently it has been found that the evolution equations describing the scale dependence of 
certain correlation functions in four- dimensional Yang-Mills theory possess a hidden symmetry. 
Remarkably enough, the emerging integrable structures are well-known in the theory of lattice 
integrable models pQ as corresponding to open Heisenberg spin magnets. In particular, the energy 
spectrum of the magnet determines the spectrum of the anomalous dimensions of the correlation 
functions in Yang-Mills theory [21 El E] • A unusual feature of these models as compared with 
conventional magnets studied thoroughly in applications to statistical physics [E] is that the 
spin operators are generators of infinite-dimensional representations of the SL(2, M) group. This 
group emerges as subgroup of the full 50(4,2) conformal group of four- dimensional Yang- Mills 
theory. 

Exact solution of the spectral problem for integrable systems with infinite-dimensional quan- 
tum space is a nontrivial task. The conventional Algebraic Bethe Ansatz (ABA) 6] is not always 
applicable to such systems and one has to rely instead on a more elaborated methods like the 
Baxter Q— operator [5] and the Separation of Variables (SoV) [7j. Being combined together, the 
two methods allow one to find the energy spectrum of the model and obtain integral representa- 
tion for the eigenstates. At present, such program has been carried out for a number of models 
with periodic boundary conditions. They include periodic Toda chain jHUH], the DST model (THj, 
noncompact closed SL{2) Heisenberg magnets [TTJIT21II3] and Calogero- Sutherland model jl4j. 
In the present paper, we apply the both methods to the quantum SL(2, K.) open Heisenberg spin 
chain. 
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A systematic approach to building quantum integrable models with nontrivial boundary con- 
ditions (including open Heisenberg spin chains) has been developed by Sklyanin ^H]- For such 
models a little progress has been made in constructing the Q— operator and the SoV representa- 
tion. One of the reasons for this is that the R— matrix formulation is more cumbersome in that 
case as compared to the models with periodic boundary conditions and, in addition, there exist 
no regular procedure for obtaining the Q— operator. 

In this paper, we construct the Baxter Q— operator and representation of the Separated vari- 
ables for the quantum SL(2, R) open spin chain. Our analysis is based on the Feynman diagram 
approach described at length in previous publications ^21 El In this approach, one realizes the 
Q— operator as an integral operator acting on the quantum space of the model and represents 
its kernel as a certain Feynman diagram. Then, various properties of the Q— operator can be 
established by making use of a few elementary diagrammatical relations. Using the obtained 
expressions, we determine the energy spectrum of the open Heisenberg spin chain in terms of the 
eigenvalues of the Q— operator and obtain integral representation for the eigenfunctions. 

The presentation is organized as follows. In Section 2 we define the open Heisenberg magnet 
with the SL(2, M) spin symmetry and review Sklyanin's formulation of the model. In Section 3 we 
construct Baxter Q— operator for the homogeneous open spin chain and establish its properties. 
In Section 4 we present an explicit construction of the unitary transformation to the Separated 
Variables for the open SL(2, M) spin chain. In particular, we calculate the integration measure 
defining the scalar product in the SoV representation and discuss its analytical properties. In 
Section 5 we demonstrate that for the open spin chain with two sites the eigenvalues of the 
Q— operator coincide with the Wilson orthogonal polynomials. Section 6 contains concluding 
remarks. Some technical details and description of the diagrammatical technique are given in 
the Appendix. 

2. Open Heisenberg spin chain 
2.1. Definition of the model 

The homogeneous open Heisenberg spin chain is a lattice model of N interacting spins S n — 
(S^, S%, S„) (with n — 1, N) described by the Hamiltonian 

N-l 

H N = Y^ H n ,n+l, #n,n+l = 2 [lp(J n ,n+l) ~ ^(2s)} , (2.1) 

n=l 

where ip(x) = d\ogT(z)/dz is the Euler ?/>— function. The pairwise Hamiltonian H n ^ n+ i defines 
the interaction between two neighboring spins S n and S n+ i. It is expressed in terms of the 
operator J n ^ n+ i related to their sum 

Jn,n+l(Jn,n+l ~ 1) = (Sn + Sn+lY . (2-2) 

The spin operators in different sites commute with each other and obey the standard commutation 
relations 

[S n i Sfc] = i£abc$nk,S n , S n = S n (s n 1) . (2-3) 

We shall assume for simplicity that the spin chain is homogeneous, s\ = ... = sn = s, with real 
s > 1/2 the same as in (|2.1j) . 
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Notice that the Hamiltonian ()2.1)1 does not involve interaction between the boundary spins 
Si and Sn- If one added the corresponding two-particle Hamiltonian H N ,i to the r.h.s. of (j2.1j) . 
the resulting Hamiltonian would define a homogeneous closed Heisenberg spin chain. The latter 
model admits solution within the R— matrix approach both by the ABA method 6} and the 
methods of the Baxter Q— operator and SoV ^3] • In the present paper we extend the analysis 
performed in the papers ^T] EI] to the case of the open spin chain and apply the method of the 
Baxter Q— operator to solve the spectral problem for the Hamiltonian ()2.1|) 

H N V q (zi, ...,z N ) = q (zi, ...,z N ). (2.4) 

Here q denotes the complete set of the quantum numbers parameterizing the energy spectrum 
and z n (with n = 1, . . . , N) are the coordinates on the quantum space V n associated with the 
nth site of the spin chain. 

The Hamiltonian Hn acts on the quantum space of the model Vn = Yln=i an d its energy 
spectrum depends on the choice of the Hilbert space V n . In what follows we shall assume that 
Vn is spanned by functions *&(zi, . . . , z N ) G holomorphic in the upper half-plane Imz n > 
and normalizable with respect to the scalar product 

(* 1 |* 2 > = Jv N z (*i(zi,... * 2 ( Zl ,...,z N ), (2.5) 

where integration measure is defined as V N z = Yln=i ^ z n with (z n = x n + iy n ) 

Vz n = d 2 z n (2 Im z n ) 2s - 2 6{lm z n ) = dx n dy n {2y n ) 2s - 2 6{y n ) (2.6) 

71 7T 

and integration in ()2.5j) goes over the upper half-plane. The spin operators S n can be realized 
on this space as differential operators 1 

Si = z 2 n d Zn +2sz w S~ = -d Zn , S° n = z n d Zn + s. (2.7) 

where = S^±iS 2 and 5° = S^. These operators are anti-hermitian with respect to the scalar 
product (1231) 

{Sn) = —Sn J (Sn)^ = ~Sn • (2-8) 

Notice that the quantum space of the model is infinite-dimensional for arbitrary finite N. For 
integer and half-integer s, the Hilbert space V n coincides with the representation space of unitary 
representation of the SL(2,"R) group of the discrete series [To] . 

2.2. R— matrix formulation 

The open SL(2,M) Heisenberg spin magnet (|2.1j) is a completely integrable model. To identify 
its integrals of motion we follow Sklyanin's approach To begin with, one defines the Lax 
operator for the SL(2,M.) magnet 

L n (u) =u + z(a.S n )=( U + f » u %o ) , (2-9) 

1 In Yang-Mills theory the spin operators H2.7[) define representation of the generators of the collinear SL(2,M) 
subgroup of the full S'0(4, 2) conformal group on the space of correlation functions (0|$ s (zin) . . . <£> s (z/vn)|0) of 
primary fields with conformal spin s and "living" on the light-cone n 2 ^ = 0. 
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where a — (01, 02, 03) are the Pauli matrices. It acts on the tensor product of the auxiliary space 
and the quantum space in the nth site, C 2 <8> V n . Taking the product of N Lax operators along 
the spin chain in the auxiliary space one defines the operator - the monodromy matrix for the 
closed spin chain 

a(u) b(u) 



Tkt(u) 



LAu) . ..L N (u) 



(2.10) 



c(u) d{u) 

which is a 2 x 2 matrix with the entries a(u), . . . , d(u) being operators acting on Vn- It satisfies 
the Yang-Baxter commutation relations 



R 12 (u - v)T N (u)T N (v) = T N (v)T N (u)Ri 2 {u - v) 



(2.1i; 



where Tn{u) = Tn(u) <8> 1 and Tn(v) — 1(g) Tn(v). The R— matrix acts on the tensor product of 
two auxiliary spaces, C 2 (g> C 2 , 

R 12 (u) = ut + iP 12 , (2.12) 

with P u being the permutation operator. The monodromy matrix for the open spin chain is 
defined as jT3j 2 



T N (u) 



T N (u) TzM 



■■n{u) J- N y—u + 1) 



■ T N (u) a 2 T{ 



N {-u) o 2 



(2.13) 



where the c- valued factor p(u) = (u—is)(u+i(s—l)) absorbs all poles of T N (u) and the superscript 
H' denotes transposition in the auxiliary space. It satisfies the fundamental "reflection" Yang- 
Baxter relation [T7| H3 d] 



2 1 

T N (v) Ri 2 {u + v — i) T N (u) Ru(u - v) 



R 12 (u - v) T N (u) R 12 (u + v-i) T N (v) (2.14) 



with the same R— matrix (J2.12)) . It proves convenient to change a normalization of IV (u) as 



T N (u) = p (-u) Tjv(-u) = T N (-u) a 2 T N (u) a 2 



It follows from (J2.13)) that Tjv(w) satisfies the relation 



A(u) B(u) 
C(u) D(u) 



T N (-u - i) T N (u) = [(u + is)(u - i(s - 1))} 2N t . 



(2.15) 



(2.16) 



The Yang-Baxter relation (|2.14|) leads to the set of fundamental relations for the operators 
A(u), . . . , D(u). For our purposes we will need only two of them 



B(u)B(v) =B(v)B(u), 

r>/ \r,S \ {U + V + i)(u-V - l) , . , 

B(u)D(v) = ; -A ; -D(v)B(u) + / 



(2.17) 



(u — v)(u + v) 

The monodromy matrix ()2.15j) satisfies the following relation 

1 



. , . u + v + i . 
A{u) H D(u) 



u — V 



goo 

u + v 



T N (u) 



2u — i 



2ua 2 T N (-u)a 2 - iT N (-u) 



(2.18) 



2 General definition of the integrable spin chain with nontrivial boundary conditions involves the boundary 
matrices K± ]J7\ El D2| • The Hamiltonian 12. 1|) corresponds to the simplest case K± = 1. 
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To verify it one starts with the definition of Tjv(w), Eq. (|2.15J) . interchanges the operators Tjv(— u) 
and T%(u) with a help of the Yang-Baxter relation ()2.11j) and uses the explicit expression ()2.12j) 
for the R— matrix. Substitution of (|2.15|) into (|2.18|) yields 

D{u) = — ^- \2uA(-u) - iD(-u)} , = B ^~ U) . (2.19) 

v ' 2u-i l v ' v n 2u + i -2u + i y J 

In the standard manner, the transfer matrix for the open spin chain t]y(u) is defined as the trace 
of the monodromy matrix (|2.15|) over the auxiliary space 

t N (u) = tr T N (u) = A(u) +D{u)=(l- D(u) + (l + D(-u) , (2.20) 

where in the last relation we took into account (j2.19|) . Following Sklyanin and making use 
of the Yang-Baxter relation (j2.14j) . one can show that the transfer matrix commutes with itself 
for different values of the spectral parameter, with the Hamiltonian (J2.1j) and with the operator 
of the total spin S = J2n=i $n 

[t N {u),tN{v)} = [t N {u),H N ] = [t N {u),S] =0. (2.21) 

The expansion of t^^u) in powers of u generates the integrals of motion of the model. One 
deduces from ()2.20j) and ()2.15|) that the transfer matrix is an even polynomial in u of degree 2N, 
t]y(—u) = t^iu), which scales at large u as t^{u) ~ 2(— l) N u 2N . In addition, it follows from 
firrSJ) that Tjy(i/2) = 11 leading to 

t N (-i/2) = 2p N (i/2) = 2{s - l/2) 2N . (2.22) 

These properties imply that t^iu) — t N (±i/2) is proportional to {u + i/2){u — if 2) 

t N (u) = (-1) N (u 2 + 1/4) [2u 2N - 2 + q 2 u 2N -' + ... + q N _ lU 2 + q N ] +2(s- 1/2) 2N . (2.23) 



Here the q— operators are given by polynomials in the spin operators S n , for instance 

q 2 = -4S 2 + 2Ns{s - 1) - X - . (2.24) 

It follows from ()2.21|) and ()2.23|) that — 1 operators ■ ■ ■ , (Jn form the family of mutually 
commuting SL(2) invariant integrals of motion. Since [Hn,S] = 0, the remaining iVth integral 
of motion is provided by one of the components of the total spin. It is convenient to choose the 
latter as iS_ = —i d Zn since its eigenvalues define the total momentum. 

Thus, the open Heisenberg spin chain is a completely integrable model and the spectral 
problem for the Hamiltonian (|2.1|) can be reformulated as the spectral problem for the transfer 
matrix 

t N (u)^ q, p (zi, ...,z N )= t N (u)y q>p (zi, ...,z N ), (2.25) 
(iS- - p)* giP (2i, . . . , z N ) = , 
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where tjv(w) is the eigenvalue of the transfer matrix f)2.20j) and q = (92, . . . , 5jv) denotes the 
eigenvalues of the integrals of motion. A general solution to (|2.25jl takes the form 

/oo 
dx e ipxo ty q ( Zl -x ,...,z N -x ), (2.26) 
-00 

where integration goes along the real axis. The eigenstate ^ q (zi, . . . , z^) has to diagonalize 
simultaneously the operators q 2 , . . . , qN- 

3. Baxter Q— operator 

To solve the spectral problem for the open Heisenberg spin chain, Eq. ()2.25|) . we apply the method 
of the Baxter Q— operator. The method relies on the existence of the operator Q(u) which acts 
on the quantum space of the model Vat, depends on the spectral parameter u and satisfies the 
following defining relations: 

• Commutativity: 

[Q(u),Q(t;)] = 0. (3.1) 

• Q - t relation: 

[Q(u),t N (u)]=0. (3.2) 

• Baxter relation: 

t N (u)Q(u) = A + (u)Q(u + i) + A_(u)Q(u-z), (3.3) 

where A±(u) are some scalar functions of u. For the homogeneous open spin chain they are given 
by 

A ± (u) = (-l) N2 ^(u±isr. (3.4) 

In this Section, we construct the operator Q(u) satisfying Eqs. (j3.1J) - ()3.3J) and discuss its prop- 
erties. 

It follows from (J3.1)) and ()3.2)1 that the Baxter Q— operator and the transfer matrix £/v(w) 
share the common set of the eigenstates 

Q(u)*,(2i, ...,z N ) = Q q (u)V q (zi, ...,z N ). (3.5) 

The eigenstates ^f q (zi, . . . , zn) are the solutions to the Schrodinger equation ()2.4|1 whereas the 
corresponding eigenvalues of the Q— operator, Q q (u), satisfy the Baxter relation ()3.3|) with the 
transfer matrix tjy(u), Eq. ()2.23|) . replaced by its eigenvalue. As we will show below, the Baxter 
Q— operator encodes information about the spectrum of the open spin chain. Namely, having 
calculated its eigenvalues Q q {u) one would be able to reconstruct the energy spectrum of the 
model E q . 
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3.1. Gauge transformations 

Our approach to constructing the Baxter Q— operator is based on the representation of Q(u) as 
an integral operator acting on the quantum space of the model 

[Q(u)*] {zi, ...,z N ) = J V N w Q u (z 1 z N \w u W N )^f(wi, . . . , w N ) , (3.6) 

with w n = w* and the integration measure defined in ([2.6)1 . To find the explicit expression for 
the kernel Q u {z\ . . . , zn\w\, ■ ■ ■ , wn) we shall explore the fact that the transfer matrix of the open 
spin chain, Eq. (|2.2U|) . is invariant under local gauge transformations of the Lax operators |H) ITT] 

L n {u) -> L n (u) = M~ l L n (u) M n+1 , (3.7) 

where M n (with n — 1, . . . , N + 1) are arbitrary 2x2 matrices with det M n ^ 0. According to 
Eqs. f|2.10|) . ()2.13|) and (|2.15|) the operators T N {u) and T N {u) are transformed under ()3.7j) as 

T N {u) -> f N {u) = M^ 1 T N {u) M N+1 , 

f N (u) -> f N (u) = Mf 1 f jv(«) Mi , (3.8) 

so that the transfer matrix t^iu) = trlV(w) stays invariant. 

The gauge rotation of the Lax operator, Eq. ()3.7j) . has been used by Pasquier and Gaudin 
to construct the Baxter operator for the periodic Toda chain jH] and it was later applied to the 
closed spin chain in Refs. fH\ I13j. In the latter case, the transfer matrix equals tiT^(u) and 
in order to preserve its invariance under the transformation (13. 8 j) one had to impose periodic 
boundary conditions M/v+i = Mi. For the open spin chain the matrices Mtv+i and Mi can be 
arbitrary, M^+i 7^ Mi. 3 In spite of this difference, many results obtained in Refs. fJJ ^] for 
the closed spin chain are applicable to the open chain. 

To begin with, let us introduce the function [TTJ ITS] 

N 

Y u (z 1 ,...,z N \w 1 ,...,w N+1 ) = f[( z k -Wk)' s ' iu {z k - w k +i)~ s+iu , (3.9) 

fc=i 

which is a (anti)holomorphic function of the complex variables z = (z\,...,zji) and w = 
(ilii, . . . , Wn+i) in the upper half-plane Im^ > and Imw n > (with w n = w*). It satis- 
fies the following relations 

b(u; w^wn+^Y^zIw) = 0, 

a(u;w 1 ,w N+ i)Y u (z\w) = (u + is) N Y u+i (z\w) , (3.10) 
d(u;wi,w N+ i)Y u (z\w) = (u - is) N Y u _i(z\w) , 

where the operators a(u) , . . . , d(u) are defined similarly to (J2.1(J|) as the entries of the gauge 
rotated transfer matrix Tn(u), Eq. ([3.8)1 . with the M— matrices given by 

mi = ( 1 v r ) . ^« = ( 1 v t +i ) • < 3 - u > 



3 Notice that the monodromy matrices Tn(u) and Tn(u) satisfy the Yang-Baxter relations, Eq H2.11(l and 
1)2.14(1 . respectively. This follows immediately from the invariance of the R— matrix 1)2.12(1 under transformations 
R^UR I/" 1 with U = M M. 
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with wn+i 7^ W\. They are given by linear combinations of the operators a(u) , . . . , d(u) , 
Eq. (|2.1()jl . with the coefficients depending on the gauge parameters wn+i and w\, which are 
identified as the right arguments of the kernel Y u (zx, . . . , Zn\wi, . . . , wjy+i)- Similar relations 

hold between the entries of the monodromy matrices Tn(u) and TV("u), Eqs. ()3.8j) and (|2.15jl . so 
that 

A(u;w 1 ) + D(u;w 1 ) = A(u) + D(u), B(it;wi) = B{u) + 0(l/w x ) , (3.12) 

where we indicated explicitly the dependence on the gauge parameter w\. Eqs. (|3.10J) play a 
crucial role in our subsequent analysis. Their derivation can be found in [TT| fTTT] . 

To proceed further let us express the entries of the monodromy ^matrix of the open chain, 
B(u) and D(u), in terms of those for the closed spin chain, a(u), . . . , d{u). One finds from ()3.8|) . 
(1TTUJ) and (ISjfej) 

B{u) = b(—u)a(u) — a(—u)b(u) , D{u) = d(—u)a(u) — c(—u)b(u) . (3.13) 



3.2. Kernel of the Q— operator 

We now turn to constructing the kernel of the Baxter Q— operator and consider the following aux- 
iliary operator G(u, v) : Vn 1 — ► Vat with the kernel given by the convolution of two Y — functions 
introduced in the subsection 13.1.1 

G UjV (z x , . . . , z N \wi, w N ) = e i7Ts( - 2N ~V jT>y2...J Vy N (3.14) 
x Y u {z\, . . . , z N \wx, y 2 , . . . , Vn-i w n ) Y v (y 2 , • • • , J/jvI^i, w 2 , ■ ■ ■ , w N ) . 

Here the integration measure T>y n is defined in (I2.fij) and the prefactor is introduced for the later 
convenience. Notice that two Y — functions in 1)3.14)1 have a different number of arguments and 
depend on the same variables W\ and wn- 

Let us demonstrate that for v = —u the operator G(w, v) satisfies the relations (j3.1j) - (|3.3D 
and, therefore, it can be identified as the Baxter Q— operator for the homogeneous open Heisen- 
berg spin chain 

Q(u) = G(u, -u) , (3.15) 



or equivalently 



Q u (z 1} . . . , z N \w u ...,w N )= e ms(1N X) {z l -w 1 ) ^(z N -w N ) 

N 



f[ I Vy n (z n _! - y n y a -{z n - y n )~ Pu {Vn ~ %-i)" a "fe - w n )^ u (3.16) 



X 

n=2 



with a u = s — iu and (3 U = s + iu. To prove ()3.15|) we apply the diagrammatical approach 
developed in Ref. (THj. In this approach, one represents the kernel G UjV (z\w), Eq. ()3.14jl . as the 
Feynman diagram shown in Fig. ^ There, the arrow with the index a that goes from y to z 
represents the factor (z — y)~ a (see Eq. (|A.1[1 ) while the black blob denotes integration over the 
position w of the corresponding vertex with the SL(2, R) measure Dw (see Eq. flA.2|) and Fig. EJ. 
The operator G(u, v) is symmetric under interchange of the spectral parameters 

G{u,v) = G(v,u), (3.17) 
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Figure 1: Diagrammatical representation of the function G U)V (z\w). For v = — u the diagram 
defines the kernel of the Baxter Q— operator, Eq. (J3.16)) . Here a x = s — ix and j3 x = s + ix. 

or equivalently G u>v (z\w) = G V;U (z\w). The proof of (|3.17j) is based on the permutation identity 
shown in Fig. |HJ Writing (3 U = (3 V + i{u — v), one replaces the left-most vertical line in the left 
diagram in Fig. ^ by two lines with the indices (3 V and i{u — v). Then, one displaces the line 
with the index i(u — v) across the diagram to the right with a help of the permutation identity 
until it merges with the right-most vertical line and changes its index to a u + i(u — v ) = a v (see 
Ref. for details). The resulting diagram coincides with the original one but with the spectral 
parameters interchanged. Furthermore, it follows from (j3.14|) and (j3.10|) that 

a(u;wi, w N ) G UtV {z\ w) = (u + is) N G u+i , u (z]w) , 

d(u)Wi,w N )G U!V (z;w) = (u - is) N G u _ itV (z;w) , (3.18) 
b(u]Wx,w N )G U)V (z-,w) = 0, 

where z — (zi, . . . , z^) and w = (wi, . . . , Wn). Here the operators a(u), . . . , d{u) depend on the 
gauge parameters vbi and wn, which coincide with the corresponding arguments of the kernel 
G UjV (zx, . . .,z N \wi, . . .,%). 

Let us demonstrate that the operator Q(u) satisfies the Baxter equation (j3.3J) . To this end, one 
examines the expression entering the l.h.s. of the Baxter equation (J3.3|) and applies Eqs. (J2.2(J|) . 
fTFJ and (ETTol) to get 

t N (u) Q(u) = (A(u) + D(u)) G(u, -u) = 

Taking into account Eqs. (|3.13|) . (|3.18|) and (|3.17|) . one finds 

D(u) G(u, —u) = d(—u)a(u)G(u, —u) — (u + is) N d(— u)G(— u, u + i) 

= (u + is) N (-u-is) N G(-u-i,u + i) (3.20) 

Substituting ()3.20|) into ()3.19|) one concludes that the operator Q(u) defined in ()3.16p verifies 
the Baxter relation (|3.3j) . In addition, one deduces from ()3.13|) and (j3.18J) that the kernel of the 
Q— operator is nullified by the operator B{u) 

B(u]w 1 )Q u (zi,...,z N \w 1 ,...,w N ) = 0. (3.21) 



2u - i ~ . . 2u + i ~ 

D(u) H z D(-u) 



2u 



2u 



G(u,-u). (3.19) 
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We will use this property in Sect. H]to construct the unitary transformation to the SoV repre- 
sentation. 

The operator Q(w) has the following properties 

• Parity: 

Q(u) = Q(-u) . (3.22) 

• Normalization: 

Q(±zs) = K. (3.23) 

• Hermiticity: 

(Q( U ))t = Q(u*) . (3.24) 

• SL(2) invar iance: 

[Q(u),£] = 0, (3.25) 

where IK is the unit operator on the Hilbert space of the model Vn (see Eq. IA.4I in Appendix IXjl 
and S = J2k=i $k is the operator of the total spin. Eq. ()3.22|) is a consequence of (j3.15J) and 
([3.17)1 . Eq. ()3.23|) follows from the fact that /3j s = a_j S = so that the corresponding lines in 
the diagram in Fig. disappear leading to drastic simplification of the kernel. Eq. ()3.24j) follows 
directly from the definition of the conjugated operator (Q(w)) . To verify (|3.25j) one notices that 
the kernel of the Q— operator, Eq. (|3.1fij) . is transformed under the SL(2,M.) transformations as 

N 

Qu(z\w') = ]J(jw k + 5) 2s ( 1Zk + 5) 2s Q u (z\w) (3.26) 

k=l 

where z' k = [az k + (3)/{^/z k + 5) and w' k = (aw k + (3)/{^w k + 5) with real a, ... ,6 such that 
a5 — /?7 = 1. 

We are now ready to demonstrate that the Q— operator (|3.15j) satisfies the relations ()3.2|) and 
(13. lj) . To verify (|3.2j) . one performs the Hermitian conjugation of the both sides of the Baxter 
equation (j3.3J) . Taking into account ()3.24|) . one finds that the r.h.s. of ()3.3|) goes into £/v(«*)Q(w*) 
whereas its l.h.s. is replaced by (tTv(-u)Q('u))^ = Q(w*)£jv(w*). Equating the two expressions one 
arrives at ()3.2j) . Finally, let us show that the operator (|3.15)l satisfies the commutativity condition 
(|3.1|) . The proof can be performed diagrammatically. To this end, one examines the Feynman 
diagram corresponding to the product Q(v)Q(u) = G(v, —v) G(u, —u) and inserts a pair of lines 
with the indices ±z(u + v) into one of the central rhombuses as shown in Fig. |21 Displacing 
the two lines horizontally in the opposite directions with a help of the permutation identity (see 
Fig.[8J) one obtains the Feynman diagrams shown in Fig.|21to the right. It differs from the original 
diagram in that various a— and (3— indices got interchanged and two additional lines with the 
indices ±i(u + v) connect the "end points", w\ with z\ and wn with z^. Taking into account 
the definition of the G— operator, Eq. ()3.14|) (see Fig. one finds that the Feynman integral 
corresponding to this diagram can be written as 



(G(u*,v*)) j G(-u,-v) (z;w), (3.27) 



where the kernel of the integral operator (G(u*, v*)y is given by (G u * )V *{wi, . . . , Wn\zi, . . . , z N ))*. 
According to (|3.17|) . the r.h.s. of ()3.27|) is invariant under interchanging u ^=> v thus proving the 
commutativity relation (|3.1|) . 
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w 1 w 2 • • • w N w 1 w 2 • • • W N 




Zl Z 2 • • • ZN Z\ Z 2 • • • Zn 



Figure 2: Diagrammatical proof of Eq. (|3.27|) . The left diagram represents the kernel of the 
operator Q(v)Q(u). The right diagram is obtained by displacing two wavy lines carrying the 
indices ±i(u + v) to the right/left with a help of the permutation identity. Here 

j3 x = a- x = s + ix. 



3.3. Contour- integral representation for the Q— operator 

In the previous subsection we constructed the Baxter Q— operator for the homogeneous open spin 
chain, Eqs. (j3.16j) . As was already mentioned, the Q— operator is diagonalized by the eigenstates 
of the model, Eq. ()3.5j) . and the corresponding eigenvalues Q q {u) satisfy (|3.3|) . 

The Baxter equation ()3.3j) is a finite-difference functional equation and its solutions are defined 
up to multiplication by an arbitrary periodic function, f(u + i) — f{u). To fix this ambiguity 
and determine eigenvalues of the Q— operator, one has to specify analytical properties of Q q {u). 
They can be identified using the following contour-integral representation for the Q— operator on 
the quantum space of the model *f?(zi,..., z?v) £ Vat 

[Q(u)V] (z h ...,z N ) = [B(s + tu,s- iu)]~ 2N+1 (3.28) 

rl N r l N 

x / n^i-^r+^-vr^ 1 / n d7 *( 1 - 7 *)* +4tt ~ lr fc" itt " 1 *(^---'^)' 

J ° n=l J( > k=2 

where B(x,y) is the Euler beta-function and the Z— coordinates are defined as 
Z t = (1- a 1 )z 1 + a 1 [r 2 z 1 + (1-t 2 )z 2 } 

Z k = (1 - a k )[TkZk-i + (1 - r k )z k ) + a k [r k+1 z k + (1 - r k+1 )z k+1 ] , (1 < k < N) , 
Z N = (1 - ct n )[t n z N -. 1 + (1 - t n )z n ] + a N z N . (3.29) 

To obtain ()3.28|) . one uses the integral representation for the Q— operator, Eq. (j3.16|) . and applies 
the identity (TO|l . 

Since the function *&(Zi, . . . , Z^) is holomorphic in the upper half-plane Imz n > 0, the 
integral in the r.h.s. of (J3.28|) is convergent inside the strip — s < Re(iu) < s in the complex 
u — plane. Analytically continuing the integral outside this strip, one finds that it contains poles 
of the order p < 2N—1 originating from integration at the vicinity of the end-points a n , t& — > 0, 1. 
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They are located at iu = ±(n — s) (with n nonnegative integer) and are compensated by the 
beta-function prefactor entering ()3.28j) . As a result, [Q(u)\E r ](zi, . . . , z^) does not have poles in 
u and, therefore, the eigenvalues of the Baxter operator, Q q (u), are entire functions of u. 

Eq. ()3.28|) also allows one to determine asymptotic behaviour of Q q (u) at large u. It is given 

by 

Q q (u) ~ u 2h [I + 0(l/u 2 )] , (3.30) 
with h nonnegative integer defining the total spin of the model 

[S 2 -(h + Ns)(h + Ns- = . (3.31) 

To establish ()3.3())1 . one verifies using 1)3.28)1 that the Baxter operator is invariant under arbitrary 
SL(2, C) transformations, in particular under the following one z i— > —i(w +i)/(w — i) 

. . . , z N ) ~ ...,%) = TTK - *T 2s * i^^ 1 ., • • • , , (3.32) 

which map the upper half-plane Imzfc > into a unit disk \wt\ < 1. The main advantage of 
dealing with functions ^{w\, . . . ,wn) holomorphic inside the unit circle is that solutions to (j3.31j) 
have a simple form in that case. Namely, the Hilbert space of the model contains the highest 
weigths which satisfy ()3.31|) and are given by homogeneous translation invariant polynomials of 
degree h, ^{w\, . . . , w^) = Ph(wi—w 2 , ■ ■ ■ , wn-i—wn). Since the Baxter operator ()3.28|) remains 
invariant under ()3.32|) . one can substitute the function &(z) in ()3.28|) by such polynomial. Then, 
^(Zi, . . . , Zjv) entering ()3.28|) becomes a polynomial in the a— and r— parameters. Integrating 
term-by-term in the r.h.s. (|3.28|) one finds that the dominant contribution at large u comes from 
terms containing a maximum number of cr's and r's. This number equals 2h and leads to the 
asymptotics ()3.30|) . 

Given that Q q {u) is an even function of u, Eq. (|3.22p . and making use of (|3.24|) . we conclude 
that the eigenvalues of the Q— operator are real polynomials in u 2 of degree h 

h 

Q q (u) = a q l[(u 2 - X 2 k ) , (Q q (u))* = Q q (u*) (3.33) 

k=l 

with the normalization constant a q fixed by the condition Q q (is) = 1, Eq. ()3.23|) . Substituting 
(|3.33p into ()3.3|) and putting u = A|, one finds that the roots satisfy the Bethe equations for 
the open spin chain |15j . 

3.4. Relation to the Hamiltonian 

Let us demonstrate that the Hamiltonian of the open spin chain, Eq. (|2.1|) . is given by a loga- 
rithmic derivative of the Baxter operator evaluated at special values of the spectral parameter 
u = ±zs. Due to ()3.23|) the expansion of the Q— operator around u = ±is can be written as 

[Q(±ts + e)V}(z) = V(z) t ie[H N mz) + 0(e 2 ) , (3.34) 

with T~Cn being some integral operator. Its explicit form can be found from the contour- integral 
representation for the Q— operator, (|3.28|) . At u — —is + e the beta-prefactor in the r.h.s. of 
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()3.28)1 vanishes as e 2 ^ 1 but it is compensated by poles coming from integration at the vicinity 
of a n = Tk = 0. Carefully separating contribution from this region, one obtains that the operator 
1~Ln entering (J3.34j) is given by the sum of two-particle integral operators 



d V "' 
H N = -i—\nQ(-is + e 
de 



e=0 

n=0 



Hn,n+1 , (3.35) 



where H n>n+ i acts on *$?(z n , z n+ i) EV n ® V n+ \ as 

f 1 dr 

[H 12 *]( Zl ,z 2 ) = \ —(l-r) 2s - 1 ^(z 12 (T),z 2 ) + ^(z 1 ,z 21 (r))-2^(z 1 ,z 2 )}, (3.36) 

Jo T 

with Zik(r) = (1 — r)Zi + TZk- It is straightforward to check that the Hamiltonian H njU+ i 
commutes with the two-particle spin S n + S n +i defined in ()2.7|) and, therefore, it only depends on 
the Casimir operator J„ )n +i, Eq. (|2.2j) . To find the explicit form of this dependence one applies 
Hi 2 to the state \I/(zi,z 2 ) = (z\ — z 2 ) h /((zi + t){z 2 + i)) h+2s with h nonnegative integer. 4 It 
diagonalizes simultaneously the Casimir operator J 12 ^ = (h + 2s) ^ and the two-particle kernel 
Hi2$ = 2[i[)(h + 2s) — ■0(2s)] 1 I r leading to the expression 

H nin+1 = 2 [V(J„,n+i) - V>(2s)] , (3.37) 

which coincides with (|2.ip . 

Eq. f!3.35|) establishes the relation between the Hamiltonian of the model ()2.1|) and the Baxter 
Q— operator. Obviously, the same relation holds between their eigenvalues 

E q = ±i—\nQJ±is + e) . (3.38) 
de e=o 

Thus, to reconstruct the energy spectrum of the model, one has to find polynomial solutions to 
the Baxter equation ()3.3j) and apply ()3.3 



4. Separation of Variables 

In this section we will construct integral representation for the eigenstates of the model, Eq. ()2.4|) . 
by going over to the representation of the Separated Variables (p, x) = (p, x%, ...,xn-i) (SoV) 

(zi, . . . ,z N ) = / d N ' 1 xfx(x) U px ( Zl ,..., z N )$ q (x). (4.1) 

Here $ q (cc) is the eigenfunction of the model in the separated variables. It is factorized into 
a product of functions depending on a single variable § q (x) ~ Q(x\) . . - Q(xn-i)- As will be 
shown in this section, Q(xk) coincides with the eigenvalue of the Baxter Q— operator. The kernel 
U PtX of the unitary operator corresponding to the SoV transformation is defined as 

U p , x (zi, . . . , z N ) = (zi, . . . , z N \p, x) . (4.2) 
4 Undcr conformal mapping (|3.32l) this state is transformed into a homogeneous polynomial of degree h, 
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We will argue below that the separated variables Xk (with k = 1, . . . , N — 1) take real positive 
values so that integration in (j4.1j) goes over x e with d N ~ 1 x = dx\ . . .dx^-i and jJ,{x) 

being a nontrivial integration measure. Eq. (j4.1)l defines the transformation $ q i— > ^^p. The 
inverse transformation looks as follows 



$ q (x) 5(p - p') = (p, x\^ q>p ) = J V N z (U p ^ x (z 1 , z N ))*^ q>p (zi, ...,z N ) 



(4.3) 



To construct the unitary transformation to the SoV representation one has to specify the 
complete set of the states \p, x) and define the corresponding kernel (J4.2)) . It is well-known 
that for the SL(2) spin chain with periodic boundary conditions, within the framework of the 
Sklyanin's approach [Zj, the basis vectors \p, x) can be defined as eigenvectors of the operator b{u) 
which is the off-diagonal matrix element of the monodromy matrix T^(u), Eq. (j2.10j) . We will 
demonstrate that the same recipe also works for the open spin chain. Namely, the basis vectors 
\p, x) in (14. 2 J) can be defined as the eigenstates of the operator B{u) entering the expression for 
monodromy matrix Tn(u), Eq. (j2.15|) . 

According to (j2.15j) . B(u) is a polynomial in u of degree 2N — 1 with operator- valued coef- 
ficients, B{u) = 2z(-l) Af - 1 S_w 2Af - 1 + . . .. In addition, it follows from (J2~TH|) that B(-i/2) = 
and, moreover, B(u)/(2u + i) is an even function of u. This suggests to remove the "kinematic" 
zero of B(u) and define the operator 

B(u) = = {-If -US- (u 2N ~ 2 + bW N -' + ... +b N ), (4.4) 

with 6 2 ,...,&7v being some (commuting) operators. Since B(u) = b(—u)a(u) — a(—u)b(u) (see 
Eq. 1CTI3I) ). one finds using (a(M)) 1 " = a(u*) and {b(u)) ] = b(u*) that (B(u)^ = -B(-u*), or 
equivalently (B(u)Y = B(u*). Thus, B(u) is hermitian operator for real u. 

Following Sklyanin we identify the eigenstates of the operator B(u) as the kernel of the 
transition operator to the SoV representation 

B(u) Up^Z!, ...,z N ) = (-l)^ 1 p (u 2 -xl)--- (u 2 - x^_ x ) U PjX (z!, ...,z N ). (4.5) 

According to ()2.17|) . [B(u), B(v)] = and, therefore, U p . x (zi, . . . , z N ) does not depend on the 
spectral parameter u. Due to ()4.4|) . the corresponding eigenvalues are real polynomials in u 2 of 
degree N — l. They can be parameterized by the total momentum p and by the set of parameters 
x = (xi, . . . ,xn-i) which are identified as the separated variables. Hermiticity of the operator 
B{u) implies that x 2 k can be either real, or can appear in complex conjugated pairs, x\ = (x 2 )*. 
We will argue in the next section that the separated variables satisfy a much stronger condition 
x\ > 0, which together with the symmetry of ()4.5|) under Xk — > — Xk allows one to assign to the 
separated variables x real positive values. This follows from the requirement that U PjX (z\, . . . , z^) 
have to be the eigenstates of the self-adjoint operator B{u) and, therefore, they have to fulfill 
the completeness condition 

r°° c ^ . 

I dpi d N ~ 1 xfi(x) (U PmSB (wi, . . . ,w N ))* U PjX {z x , . . . ,z N ) = T\^{zn\w n ) (4.6) 

Jo ^ 

where lK(z|w) = e ms (z — w)~ 2s is the kernel of the identity operator (see Eq. (|A.4jl ). 
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Figure 3: The diagrammatical representation of the function A x (zi, . . . , zn\w 2 , ■ ■ ■ , iujv). 

The diagonal element D(±x k ) of the monodromy matrix (|2.15|) acts on U PjX (wi, . . . ,wjy) as 
a shift operator 

D(±x k )U PjX (zi, ...,z N )= 5(±Xk)Up t B±i ek (zi, ...,z N ). (4.7) 

Indeed, taking v = ±x k in the second fundamental relation (j2.17|) and applying its both sides 
to U PiX (zi, . . . ,Zn), one arrives at ([4.7)1 . The scalar factor S(xk) depends on the normalization 
of U P)SB (zi, . . . , zn)- Applying U P)X to the both sides of (|2.16j) and taking u = x k one finds that 
S(x k ) satisfies the relation 

5(x k )5(-x k -i) = [(x k + is)(x k + i(l - s))] 2N . (4.8) 

In (14. 7|) it was tacitly assumed that the function U p>as (zi, . . . ,zjf) can be continued to complex 
x. Notice that U PtX ± iek (zi, . . . , z N ) is not the eigenfunction of the operator B(u) even though it 
satisfies the differential equation (|4.5j) . 

4.1. Transition kernel 

Solving the spectral problem (|4.5|) . we follow the approach developed in Ref. jT3] in application 
to the closed spin chain. To begin with, we notice that the differential equation ()4.5|) is equivalent 
to the system of iV equations 

iS-Up^z) = P U p>x {z), B(±x k )U p , x (z) = 0, (fc = l,..JV-l). (4.9) 

Let us consider the second relation and compare it with a similar relation (j3.21|) for u = ±x k . 
Sending the gauge parameter w% in (|3.21|) to infinity and taking into account ()3.12|) one finds 
that B(±x k ) annihilates the following function 

A x (z 1 ,...,z N \w 2 ...,w N ) = lim wl s Q x (z x ,...,z N \wi,...,w N ). (4.10) 

MIX— >0O 

Here the additional prefactor is inserted to make the limit finite 

A u (zi, . . . , z N \w 2 ■ ■ ■ ,wn) = e^ 2 ^ 1 ) J Vy 2 ...Vy N {y N -w N )-^{z N -w N )- a - (4.11) 

N N-l 

x n^^- 1 ~ yk)~ au ( z k - Vk)' Pu n (v* - w n y pu {y n +i - w n y au . 

k=2 n=2 
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As before, it is convenient to represent this expression as the Feynman diagram shown in Fig. El 
It differs from the Feynman diagram for the Q— operator (see Fig. in that two lines attached 
to the vertex w\ are removed. 

By the construction, the A— function satisfies the relation 

B(±x k )A Xk (z 1 ,...,z N \w 2 ...,w N ) = 0. (4.12) 

Let us introduce the integral operator An{x) with the kernel given by (|4.11|) . It maps a function 
of N — 1 variables ^n-i(w 2 , • • • , wn) into a function of N variables ^n{zi, • • • , z N ) 

■V N ( Zl ,...,z N ) = [A N (u)V N - 1 ](z 1 ,...,z N ) (4.13) 

Vw 2 ... / Vw n A u (z 1 ,... j z n \w 2 ...,wn)'^n-i(w2,...,w n ), 



The operator Ajy(x) defined in this way has a number of remarkable properties: 
• Parity: 



Commut at i vity : 



Baxter relation: 



A N (x) = A N (-x) (4.14) 

A N (X 1 )A N - 1 (X 2 ) = A N (x 2 )A N - 1 (x 1 ) (4.15) 



t N (x) Ajv(ar) = A+(x) A N (x + i) + A_(x) A N (x - i) (4.16) 

• Exchange relation: 

A j N (x)A N (y) = tp(x, y) ■ A JV _ 1 (y)A] v _ 1 (x) , (4.17) 
where x ^ y and the scalar function cp(x, y) = (p(y, x) is defined as 

ip(x, y) = e 4i7rs a(a x , a y ) a((3 x , (3 y ) a{/3 x , a y ) a(a x , (3 y ) (4.18) 
with a x = s — ix and /3 X = s + ix, 

a(a ,,) = 9 -~. r(a+ r f-^ r(2S > . (4,9) 

The following comments are in order. 

Eq. 1)4.14)1 follows from the parity property of the Baxter Q— operator, Eq. (J3.22)) . The proof 
of (j4.15j) can be performed diagrammatically with a help of the permutation identities (see Figs. |H1 
andEJ). It goes along the same lines as the proof of commutativity property for the Q— operator 
presented at the end of Sect. 3.2. Eq. ()4.16|) follows immediately from (|4.10|) and ()3.3|) . The proof 
of the exchange relation (|4.17|) is illustrated in Fig. HJ The product Ajy(a;)Ajv(z/) corresponds to 
the left diagram in Fig. 0] The left-most vertex in this diagram can be integrated out with a help 
of the chain relation (see Fig. Q) producing a single line with the index a x + (3 y — 2s = —i(x — y). 
Then, one moves this line horizontally to the right of the diagram by applying the permutation 
identity (see Fig. |HJ). Repeating the same steps for the resulting diagram one finally arrives at 
the right diagram in Fig. 0] with the additional prefactor f)4. 18j) . 
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Figure 4: Exchange relation. Here and P x = s + zx. 

Taking into account the properties of the A— operator, it becomes straightforward to write a 
general solution to the system (|4.9j) 

U p>x (z u ...,z N )= p Ns ^ 2 J Vw N e™ U x (z; w N ) , (4.20) 

where U x (z; wn) is factorized into the product of iV — 1 operators 

U x (z; w N ) = [A N (xi) AiY_i(x 2 ) . . . A 2 (xv-i)] {z u . . . , zjvI^at) , (4.21) 

with z — (z\, . . . , zn) and the additional factor p^- 1 / 2 introduced in ()4.20j) for the later conve- 
nience. Indeed, the first relation in (|4.9jl is satisfied due to invariance of (j4.21|) under translations 
z k - z k + e and — > it; at + e with e real. It follows from (|4.14j) and (|4.15j) that U x (z; wn) is an 
even symmetric function of Xi, . . . ,Xtv-i- Since B(±xi) U x (z;w N ) = by virtue of (|4.21|) and 
([4. 12)1 . the second relation in (j4.9J) is fulfilled for arbitrary k. Notice that the kernel U x (z] wn) 
satisfies a multi-dimensional Baxter relation 

t N (x k ) U x (z; w N ) = A + (x k )U x+iek (z; w N ) + A_(x fc )E7 (B _ feh (z; w N ) , (4.22) 

where denotes a unit basis vector in the x— space, x = Yl,k x k e k- Eq. ()4.22|) follows from the 
similar property of the A— operator, Eq. (|4.16p . and the symmetry of the kernel under permuta- 
tions of x— variables. 

Eqs. (|4.2U|) and (|4.21|) define the transition kernel U PtX (zi, . . . , z N ) to the SoV representation 
for the homogeneous open spin chain. Remarkably enough, these expressions have the same 
form as for the closed SL(2) spin chain ^H]- The only difference between the two cases is in the 
definition of the A— operator. Diagrammatical representation for the transition kernel (|4.21j) is 
shown in Fig. The corresponding Feynman diagram has a pyramidal form which reflects the 
structure of the kernel (J4.21|) . It consists of (N — 1)— rows with each row representing a single 
A— operator. 

It remains to verify that the kernel ()4.20|) satisfies for real x the completeness condition (|4.6|) . 
As we will show in Sect. 5, the transformation to the SoV representation for N = 2 open spin 



17 




Z2 



Z 3 



z 4 



Figure 5: The diagrammatic representation of the kernel Ux{z,wn) for N = 4. 

chain coincides with the Fourier- Jacobi transform (see, e.g. Ref. JE]). Then, reality condition 
for the separated variable x and completeness condition for U p x (zi, z 2 ) follow immediately from 
the properties of the Fourier- Jacobi transform. For N > 3 some arguments will be presented in 
the next subsection. 



4.2. Integration measure 

Let us demonstrate that the transition kernel defined in Eqs. f|4.20|) and (j4.21|) satisfies the 
orthogonality condition 



(p', x'\p, x) 



V zU PiX (zx, ... } z N )(U p/>x >(zi, ...,z N )) 



S(p-p') {5(x — x') H } 



(N-iy. 



(4.23) 



and calculate the integration measure n(x). Here 8{x — x') = Ylk=i $( x k ~ x 'k) an( ^ x = 
(x\, . . . , xn-i) take positive real values, Xk > 0. Ellipses denote the terms with all possible 
permutations inside the set x = (xi, . . . , xzv-i)- 

The calculation of (|4.23|) repeats similar analysis for the closed spin chain described at length 
in Ref. [T3]. Substitution of flOOD into (|Q3J) yields 



(p',x'\p,x) = (pp') Ns - 1/2 J Vw N e ipWN J Vw'j 



, N e- i * iB '»(w' N ,x'\w N ,x) 



(4.24) 
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Figure 6: The scalar product (w',x'\w,x) for N = 2 

where the notation was introduced for the ket-vector (z 1 , . . . , z n \wn, x) = U x (z; wn), or equiv- 
alent ly 

\w N , x) = A. N (x 1 )A N ^ 1 (x 2 ) . . . A 2 (x N _ 1 )\w N ) , (4.25) 

with \wn) being a "single-particle state". We recall that the operator A k (u) maps (A; — 1)— particle 
state into A;— particle one, so that a composition of the A— operators in ()4.25|) produces the 
N— particle state. Calculating the scalar product (w' N , x'\wn, x) one applies systematically the 
exchange relation (|4.17jl and obtains 

(w' N ,x'\w N ,x) = c(x,x')(w' n \(aI(x' n _ 1 )A 2 (x 1 )\ . . . (^Al(x' 1 )A 2 (x N _ 1 )^\w N ) , (4.26) 

where c(x,x') = r|i<j,fc<jv-2 ip(xj,x' k ). Notice that the exchange relation ()4.17|) holds only for 

j~\-k<N— 1 

x 7^ y. Therefore, calculating ()4.26|) we have tacitly assumed that Xj ^ x' k for j + k < N — 1, or 
equivalently that all factors (p(xj,x' k ) are finite. For > 3 the matrix element entering ()4.26|) 
can be represented as follows 

J T>Wx... j T>W N ^2(W' N ,X' N _ 1 \W 1 ,X 1 )(W 1 ,X' N „ 2 \W2,X 2 ) ■ ■ ■ (wAr_ 2 ,xi|wAr,XAr_i) , (4.27) 

where (w',x'\w,x) = [A|(x') A 2 (a;)](w / ; w). Thus, the calculation of the scalar product (|4.24jl for 
arbitrary A^ is reduced to the calculation of (w',x'\w,x) at N = 2. Given that the separated 
variables at A^ = 2 take real positive values we deduce from (|4.26jl that the same holds true for 
arbitrary N. 

To calculate the scalar product at A^ = 2 we apply the diagrammatical approach of Ref. [TB*] 
and represent the matrix element (w',x'\w,x) as the Feynman diagram shown in Fig. El One 
expects from (|4.2Hjl that (w', x'\w, x) ~ 8(x— x'), so that the scalar product (w', x'\w, x) should be 
understood as a distribution. To find its explicit form we regularize the corresponding Feynman 
integral by introducing a small parameter e and shifting the indices of two lines as indicated in 
Fig. |U1 Under such regularization, the Feynman integral remains finite at x = x' and it can be 
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calculated exactly with a help of the chain relation and permutation identity (see Figs. and |EJ) • 
The calculation is straightforward and some details can be found in Ref. ^H|- Going over to the 
momentum representation and taking the limit e — > 0, one finds (for x , x' > 0) 

/ Vw e ipw (w', x'\w, x) = 2tt e ipw ' p~ 2s T 5 (2s) lEMg + g))! S ( x _ x ') . ( 4 . 28 ) 

J v ; |r(s + 2x)r(s + 2x')| 4 

Substituting (|4.2fij) and (|4.27jl into (|4.24jl and taking into account (|4.28jl one obtains after some 
algebra 

N-l N-l _ . . N.-I-2JV 

(p',x>\p,x) = (2nf-^(p-p')l[5(x k -x' N _ k ).T N (2s)ll — 



k=l k=l 



T(2s) 



x 



-l 

f] ^^sinh 7i(x k + Xj) ] [ ^-^sinh 7t(xa; — Xj) I . (4.29) 

^<j<k<N-l 71 l<j<k<N-l 71 



We recall that the calculation was performed under assumption that Xj ^ x' k for j + k < 
N — 1. Since the kernel U PiX is a symmetric function of x, (p',x'\p,x) should possess the 
same property. This allows one to relax the above assumption and replace the product of delta- 
functions rife^i 1 &{ x k— x ' N _ k ) i n the r.h.s. of ()4.29|) by the sum ^ 5 S(x-Sx') over all permutations 
inside the set x' = (V 1? . . . , x^ r _ 1 J . 

Matching (J4.29|) into (|4.23|) . one finds the expression for the integration measure in the SoV 
representation 



~ (iV-l)!(27r)^- 1 11 



T(s — ix k )T(s + ixk) 1 



(4.30) 



x 2 - x 2 2x k 
x jQ — 2 2 ^ [cosh(27TXfc) — cosh(27TXj)] Y\ sinh(27rxA;) . 

l<j<k<N-l k=l 

This expression has the following properties. As expected, fi(x) is an even function of the 
separated variables. It takes nonnegative values for real x = (xi, . . . , xn-i) and vanishes on the 
hyperplanes Xj = x k . 

After analytical continuation to complex x, the measure n(x) becomes a meromorphic func- 
tion of x k (k — 1, . . . , N — 1) with poles of the order 2N located along the imaginary axis at 
Xk = ±z(s + n) with n G N. The measure decreases exponentially fast when one of the separated 
variables, say x k , goes to infinity along the real axis 

H(x) ~ e- 2 ^! xl Ns ~ 3 . (4.31) 

One verifies that the measure (|4.3U|) satisfies the functional relation 

ji(x + ie k ) x k + i ( Xk + is \ 2 y-r Xk — Xj + % x k + Xj + i 



fi(x) x k \x k + i(l-s)J L J- x k -Xj Xk + xj 



n Xk Xj -|- Xk -I- Xj T I s 
nr. T" • T*i —I— 1" - 



with e k defined in (|P2j) . 
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It is instructive to compare 1)4.30)1 with a similar expression for the integration measure for 
the closed spin chain [T3j 

JV-l N-l 

Hd(x) = ] [ (x k - Xj) sinh(7r(x fe - Xj)) ] [ [T(s + ix k )T(s - ix k )] N . (4.33) 

i.fc=i fc=l 

One observes that /x c i(x) enters as a factor into the expression for fi(x), Eq. ()4.30j) . 
4.3. Eigenfunctions in the SoV representation 

The eigenfunctions ^ q>p (zi, . . . ,z N ) are orthogonal to each other for different sets of quantum 
numbers with respect to the SL(2) scalar product ([2.5)1 . In the SoV representation the same 
condition looks as follows 

(*«VI* W > = ($ q '\$ q ) So J(p-p') = 5(p-p f )5 q , q ', (4.34) 
where the scalar product in the SoV representation is given by 

(3y|$ 9 ) SoV = / d N - 1 xn(x)(<$> q ,(x 1 ,...,x N - 1 ))*<$> q (x h ...,x N _i) . (4.35) 

We recall that the momentum p takes real positive values whereas the spectrum of the integrals 
of motion q = (q 2 , . . . , <7jv) is discrete |£]. 

To define the eigenfunction in the SoV representation, <& q (x), we substitute \E f q iP (zi, . . . , z^) 
in (j2.25j) by its integral representation (|4.1|) . Following the standard procedure [7j and making 
use of Eqs. ()4.22)) and (|4.32|) one can show that $ 9 (a;) satisfies the (N — 1)— dimensional Baxter 
equation 

t N (x k ) ® q (x) = A + (x k )$ q (x + ie k ) + A^(x k )$ q (x - ie k ) , (4.36) 

where tN{x k ) is the eigenvalue of the transfer matrix, Eq. ()2.25jl . As before, to solve this equation 
one has to specify additional conditions for Q q (x). 

Using (j4.3j) one can show that & q {x) is a polynomial in x = (xi, . . . , sc/v-i). The proof goes 
along the same lines as analysis of analytical properties of the Baxter Q— operator in Sect. 3.3. 
Namely, substituting the expression for the kernel (J4.2U)) into ()4.3j) and applying (|A.5|) . one can 
express the r.h.s. of ()4.3|) as a nested contour integral. Analytical properties of the function 
are in the one-to-one correspondence with the properties of this integral. 

It is easy to see that polynomial solutions to ()4.36J1 can be represented in the factorized form 

$ q (x) = c q Q q (x 1 )...Q q (x N - 1 ), (4.37) 

where Q q (x) is the eigenvalue of the Baxter Q— operator, Eq. ()3.33)) . and the coefficient c q is fixed 
by the normalization condition (|4.34j) . Substituting 1)4.37)1 into ()4.35)1 and taking into account 
that Q q (x) is a real function, Eq. 1)3.33)1 . we find that the solutions to the Baxter equation satisfy 
the orthogonality condition 

» AT-l 

/ d N - x xii{x)X\Q q ,{x k )Q q {x k ) ~5 q j, (4.38) 
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with the measure given by 1)4.30)1 . 

Thus, having determined the eigenvalues of Baxter operator Q q (x) one would be able to 
restore both the energy and the corresponding eigenf unction, Eqs. 1)3.38)1 and 1)4.1)) . respectively. 
The solutions to the Baxter equation for the SL(2, R) open spin chain have been studied in jl]. 
It turns out that the ground state Sl p (z) of the model can be found exactly. This state has the 
total SL{2) spin h = 0, Eq. 1)3.31)1 . and has the form 5 

2s-l (■ N 



k=l 



Indeed, fl p (z) diagonalizes simultaneously the operators of two-particle spins, ( J n , n +i~ 2s)Q p (z) = 
0, Eq. (J2.2J) . and leads the energy E q = 0. It is interesting to notice that Q p (z) is also the ground 
state of the homogeneous SL(2,M) closed spin chain |13j . 

The state ()4.39)) diagonalizes the Baxter operator Q(u)\Q p ) = |fi p ), or equivalently Q q (u) = 1. 
As a consequence, it admits the following integral representation (see Ref. ^3j) 6 

QJz) = p Ns ~ 1 / 2 e -M2AT-i) f d N-i x , x s v _ ( 4 40 ) 

Let us calculate the SL{2) scalar product (Q P \Q P >) and use two different expressions for fl p (z), 
Eqs. ()4.39)) and ()4.40|) . Equating the two expressions, one finds that the integration measure 
satisfies the normalization condition 

dxu(x) = — - — -. (4.41) 
-i ^ V ; T(2iVs) V ; 

The relation ()4.40)) allows one to establish the equivalence between the SoV and the ABA 
methods for the SL(2, E) open spin chain. Let Ai, . . . , be the Bethe roots, or equivalently 
zeros of the polynomial Q q (x), Eq. ()3.33)) . Applying B(Xi) . . . B(Xh) to the both sides of (j4.40j) 
and taking into account ()4.5j) we obtain 

» N-l 

t> p>q (?) = B(X 1 ) . . . B(X h ) n p (z) = c(p) / dx fi(x) U p , x (z) TT Q q (x k ) , (4.42) 



where c(p) is the normalization factor. In ()4.42)) . the first relation coincides with the ABA 
representation for the eigenstate of the model while the second one defines the same eigenstate 
in the SoV representation. 

The explicit form of the eigenfunctions in the SoV representation (J4.37)) suggests that there 
exists a relation between the Baxter Q— operator and the transition kernel U PjX (z) [2D]- in the 
case of the closed spin chain it has been established in Ref. ^3]- It turns out that this relation 
is universal and it also holds for the open spin chain 

U p , x (z) = [Q(xt) ...Q(x N - 1 )e p )(z 1 ,...,z N ). (4.43) 



5 In the terminology of the Algebraic Bethe Ansatz, fl p (z) is a pseudovacuum state. 

6 Notice that this relation takes the same form both for the open and closed spins chain whereas the expressions 
for the integration measure and the transition kernel to the SoV representation are different in the two cases. 
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Here Q p (z) is a certain a;— independent function of z, which does not belong to the quantum 
space of the model. 7 Since Q(is) = K, the function Q p (z) is equal to U PjX (z) for special values 
of the x— variables, X\ = . . . = Xn-i = is that we denote as U p ^ s (z). The expression for U P)is (z) 
can be easily obtained from diagrammatic representation of the kernel (see Fig. EJ) 

%{Z) = U p , is (z) = pNs-l/2 e ins{N-l) jpz N _ (4 _ 44) 

Notice that Q p {z) depends only on a single variable z^ and, therefore, it is not normalizable with 
respect to the SL{2) scalar product (|2.5I) . 

The proof of (|4.43|) can be performed diagrammatically and it repeats similar analysis in 
Ref. [13 . Another way to verify (j4.43j) is to use the following identities 

J V N wQ u {z\w) e tpWN = e~ t7TS J Vw n K(zn-u z N \w n ) e tpWN (4.45) 
J V N w Q u {z\w)A v (w k , ■ ■ .,w N \y k+1 , ...,y N )= e~ l7TS [A k+1 (u)A k (v)](z k -i, ...,z N ; y k +i, --^Vn), 

with z — (zi, . . . , Zn), w = (wi, . . . , wn) and V N w = Yln=i ^ w n- We recall that the A— operator 
was defined in Eqs. (j4.13|) and ([4. II)) . To derive ()4.45j) one substitutes the expression for Q u (z\w), 
Eq. (j3.16|) . and integrates over "free" vertices w\, . . . ,w k -i with a help of the identity (|A.5|) for 
ty(w) = 1. Eqs. (j4.45|) can be rewritten symbolically as 

Q(u)\e p ) = e-™ A 2 (u)\Q p ) , Q(u)A k (v) = e~™ s A k+1 (u)A k (v) . (4.46) 

Applying these relations one verifies that (|4.43|) coincides with ()4.20|) and (j4.21|) . 



5. Relation to the Wilson polynomials 



In this section, we consider the open spin chain with N = 2 sites. We will demonstrate that in 
that case the eigenvalues of the Baxter Q— operator are given by the Wilson polynomials [21] 
and the unitary transformation to the SoV representation coincides with the Fourier- Jacobi 
transformation [TT?] . 

At N = 2 the Hamiltonian of the open spin chain (j2.1j) equals Ti.2 = H\i = 2[^( J12) — ip(2s)]. 
Its eigenstates are uniquely fixed by the values of the momentum p and the total SL(2) spin h, 
Eq. ()3.31|1 and are given by 



*p,a(zi,z 2 ) 



p 



2s-l 



T(2s) 



Vw e ipw 



zi - w) 2s+h (z l - w) 2s + h ' 



(5.1) 



Substituting ()5.1|) into ([3.28)1 . one finds the eigenvalues of the N = 2 Baxter operator after some 
algebra as 8 



Qh{u) 



4-^3 



-N,N + 4s-l,s + iu,s 
2s, 2s, 2s 



iu 





r T{2s) 1 


0- 


J(2s + A^)_ 



W N (u 2 , s, s, s, s) 



(5.2) 



7 In Ref. ^21; the function O p (z) was defined as a limiting case of the state \flw ,w N ) belonging to the Hilbcrt 
space of the model. 

8 As was explained in Sect. 3.3, due to the SL(2) invariance of the Baxter operator, one can calculate Qh(u) 
by substituting ^(zi, Z2) = {z\ — z-^f 1 into 
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where Wn(u 2 ) is the Wilson polynomial 21 . It is interesting to note that the solution of the 
N = 2 Baxter equation for the closed spin chain is given by the continuous Hahn polynomials 
(see e.g. Ref. ^H])- The polynomials Qh{x) are orthogonal on the half-axis x > with respect 
to the scalar product ()4.38|) with the measure ()4.30|) given by 



Hn=2{x) = — 



T 4 (s + ix) 



T(2ix)T 3 (2s) 



(5.3) 



This property is in a perfect agreement with the orthogonality condition for the Wilson polyno- 
mials |2"T] . 

The eigenvalue of the N = 2 Hamiltonian corresponding to (J5.1)) can be calculated either by 
replacing the operator J\ 2 in the expression for H 2 by its eigenvalue (Ji 2 — h — 2s)^ Pt h = 0, or 
by applying (|3.38|) . In this way, one obtains 

E h = ±i^- In Q h (±is + e) = 2 [^{h + 2s) - ^{2s)\ . (5.4) 

u€ e=0 

The ground state corresponds to h — 0. 

Let us examine the unitary transformation to the SoV representation (j4.1j) for N = 2. It is 
defined by the transition kernel, Eqs. (|4.2()jl and 1)4.21)1 . which looks at N = 2 like 



Up^zu z 2 ) = p 2s -^ 2 J Vw 2 A x {z u z 2 \w 2 ) , (5.5) 
where A x (zi, z 2 \w 2 ) are given by 1)4.11)) 

A x (z 1 ,z 2 \w 2 ) = e 3i ™ [ Vy 2 (y 2 -w 2 )-P*(z 2 -w 2 )- a *(z 1 -y 2 y a *(z 2 -y 2 )-P% (5.6) 



with f3 x = s + ix and It is convenient to transform U p ^ x (zi, z 2 ) to the momentum 

representation. 

For arbitrary function ty(zi,z 2 ) G V 2 this transformation is defined as 

*(z h z 2 ) = — — / dp 1 dp 2 e i ^ + ^(p lP2 ) s - 1 / 2 y( PhP2 ), (5.7) 

r(2s) j 

where the additional factor (pip 2 ) s ~ l l 2 was introduced to simplify the expression for the scalar 
product ()2.5)) in the momentum representation (see Eq. (|A.3|) ) 



(*i|* 2 > = J d Pl dp 2 (*i(pi,Pa)y* a (Pi,P2) • (5-8) 

In particular, the N = 2 eigenstates (j5.1j) are given in the momentum representation by the 
Jacobi polynomials 9 

*M(pi.ftO = ofc*(p-Pi -P2)(pm) s - 1/2 (pi + P 2) h v? s -^-v(Pl^Pi) , (5.9) 

\Pi + Pi J 



9 This expression is well-known in QCD as defining the conformal operators built from two fields with the 
conformal spin s. 
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where ah = i~ h ~ 4s h\T(2s)/T 2 (h + 2s) and delta-function ensures the momentum conservation. 
Applying (|A.1J) and performing integration in (|5.5jl . one finds that in the momentum represen- 
tation the N = 2 transition kernel is given by 



=Hp-p 1 -p 2 )^^(^) 1 ' 2 ( p A\fJ s -' x ; 3 + - -»). (5.io) 

\T(s + tx)\ 2 \P1P2 J \PiJ V 2s Pi J 

It defines the SoV transformation ^/ p (pi,p 2 ) i— > 

§(x)8{p-p') = dp 1 dp 2 (u p%x (p 1 ,p 2 j S j * p (pi,P2), (5.H) 

with p,x > 0. Substituting (|5.1()|) into this relation and introducing notations for ^(pi,^) = 
5(p — px — p2)p -1 ' 2 £ s-1 / 2 (l + 0/(0 w hh £ = p 2 /pi, one finds that (|5.1H1 is reduced to 

*(*) = e*" r dZe s -\FJ S - iX ; SJriX (5-12) 

This relation defines the map /(0 l— > "^K^), which is known as the Fourier- Jacobi or the index 
hypergeometric transform [T9"] . Then, the unitarity of the SoV transformation at N = 2 follows 
from the similar property of the transformation (J5.12j) . 

Let us replace ^ p (pi,p 2 ) in (|5.11|) by the N = 2 eigenstate (|5.9j) . According to (|4.37|) and (|4.3|) . 
the corresponding eigenfunction in the separated variables is given by the Wilson polynomial 
= Qh(x) ~ Wft(x 2 , s, s, s, s). Then, f!5.12|) leads to a known representation for the Wilson 
polynomials as the index hypergeometric transform of the Jacobi polynomials [Hi] . 



6. Concluding remarks 

In this paper we have constructed the Baxter Q— operator and the representation of the Sep- 
arated Variables for the open homogeneous SL(2, R) spin magnet. Our analysis relied on the 
diagrammatical approach developed in Refs. ^21 i n application to the closed spin chain. In 
this approach, one represents the kernels of the relevant integral operators (Q— operator, the 
transition kernel to the SoV representation) as Feynman diagrams and establishes their various 
properties with a help of a few simple diagrammatical identities. 

We found that the Feynman diagrams for the Q— operator and the transition kernel to the SoV 
representation have a remarkably simple form (see Figs. Q and EJ). In the latter case, the diagram 
reveals a universal pyramid-like structure which has been already observed for various quantum 
integrable models like periodic Toda chain [9,, closed SX(2,IR) and SX(2,C) spin chains [T2l ITB] 
and Calogero-Sutherland model ^3]. This structure is a manifestation of a general factorization 
property (|4.21j) of the transition kernel to the SoV representation. Namely, the kernel is factor- 
ized into the product of A— operators each depending on a single separated variable. The only 
difference between the models mentioned above resides in the explicit form of the A— operator. 
The latter can be obtained as a certain limit of the Q— operator leading to the expression for the 
transition kernel to the SoV representation as the product of the Q— operators projected onto 
a special reference state. Another advantage of the diagrammatical approach is that it offers a 
simple regular way for calculating the integration measure in the SoV representation. 
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We found that there exists an intrinsic relation between the open spin chains and Wilson 
polynomials j2I] • The latter occupy the top level in the Askey scheme of hypergeometric orthog- 
onal polynomials (221- These polynomials define the eigenvalues of the Baxter operator for open 
spin chain with N = 2 sites [3]. 

It is straightforward to extend our analysis to the case of inhomogeneous open SL(2, M.) spin 
chains. One can show that for such models the transition kernel to the SoV representation is 
given by the same pyramid-like diagram shown in Fig. El with the only difference that both the 
indices attached to various lines and the integration measure corresponding to internal vertices 
should be modified appropriately. Another interesting possibility could be to consider an open 
spin chain with the SL(2, C) symmetry. Such models naturally appear in high-energy QCD as 
describing Regge singularities of scattering amplitudes with meson quantum numbers (22]. In 
that case, the quantum space of the model does not possess the highest weight (pseudovacuum 
state) and, therefore, the Algebraic Bethe Ansatz is not applicable. 
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A Appendix: Feynman diagram technique 

Here we collect some useful formulae for the SL(2, M.) integrals. Their derivation can be found 
in Refs. [□]. 



Propagator: 

i Q — ina/2 roc 

a 1 e ' / 



w z (z-w) a r(or) .y n 

'Chain relation" (see Fig. Ej): 



dpe ip{z -^ p a -\ (A.l) 



J Vw{z-w)- a {w -v)-P = a(a,p)(z-v)- a - l3+2s , (a + f3 ^ 2s) (A.2) 



Delta-function relation: 



Identity operator: 



J Vw e ipw - ip,iD = 5{p - p') p 1 - 23 ■ r(2s) , (A.3) 



[K = J Vw ^ S _ )2 ^ {w) = ^{z). (A.4) 

Contour-integral representation: 
[Vw = T ^ /V r q*-i (1 _ r)fe -i^ 



(rzi + (1 - t)z 2 ) . (A.5) 
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Figure 7: Chain relation. 



Fourier integral: 

[z — w) a r(aj 
• Permutation identity (see Figs. |S1 and EJ): 

(z 2 - v 2 y^I(z, v;x,y) = {z x - vtf^Hz, v;y,x), (A.7) 

where z = (zi,z 2 ), v = (v 1 ,v 2 ) and 

I(z,v;x,y) — I Vw-, — — -, . , „ . rr — • (A. 8) 

J (w — v i) ax (w — v 2 ) Px [z x — w) p v [z 2 — w) a y 

In these relations, and f3 x = s + ix for arbitrary x, the a— function is defined in (j4.19|) 

and the integration measure Vw is given by ()2.6j) . w = w* and p > 0. 




Figure 8: Permutation identity. 




Figure 9: Special case of the permutation identity. It is obtained from Figure |H1 by sending one 
of the external points to infinity. 
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